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Model-Independent Properties of the B-Meson Distribution Amplitude 
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The operator product expansion is used to obtain model-independent predictions for the first 
two moments of the renormalized _B-meson light-cone distribution amplitude cf>+(uJ, (jl), defined with 
a cutoff ui < Auv- The leading hadronic power corrections are given in terms of the parameter 
A = mj- nib- From the cutoff dependence of the zeroth moment an analytical expression for the 
asymptotic behavior of the distribution amplitude is derived, which exhibits a negative radiation tail 
for uj ^> fi. By solving the evolution equation for the distribution amplitude, an integral representa- 
tion for 4>+(lo, h) is obtained in terms an initial function ^(lo, fig) defined at a lower renormalization 
scale. A realistic model of the _B-meson light-cone distribution amplitude is proposed, which satisfies 
the moment relations and has the correct asymptotic behavior. This model provides an estimate 
for the first inverse moment and the associated parameter As. 

PACS numbers: 12.38.Cy,12.39.Hg,12.39.St,13.25.Hw 



I. INTRODUCTION 



Exclusive decays of B mesons such as B — > irlv and 
B — > 7T7T, ttK are important tools to search for physics 
beyond the Standard Model as well as to measure fun- 
damental parameters in the flavor sector. In processes 
where large momentum is transferred to the soft specta- 
tor quark via hard gluon exchange, the i?-meson light- 
cone distribution amplitude (LCDA) enters in the pa- 
rameterization of hadronic matrix elements of bilocal cur- 
rent operators 0. The past few years have seen a lot of 
progress in the theoretical framework for the analysis of 
exclusive B-meson decays, mainly based on QCD factor- 
ization theorems 

HHHK and perturbative QCD meth- 
ods 

HUSH- However, in many cases the extraction of 
important physics from experimental data is still limited 
by theoretical uncertainties, often due to our ignorance 
of the functional form of the S-meson LCDA and other 
hadronic matrix elements. For exam ple, using the soft- 
collinear effective theory 0, ITTl El Mill, the large- 
recoil heavy-to-light form factors relevant to weak B de- 
cays have been studied at leading order in a l/E expan- 
sion |l5lll6lll7| . The analysis of spin-symmetry violating 
contributions to these form factors, in particular, relies 
on knowledge about the S-meson LCDA ElEHMlHI. 

In spite of the importance of the B- meson LCDA, so far 
most studies of its properties have been limited to model- 
dependent analyses based on QCD sum rules |ll l'2'2ll'23| . 
In the present work, we employ the operator product ex- 
pansion (OPE) to explore some model-independent prop- 
erties of the LCDA. We calculate the first two moments of 
the distribution amplitude, derive its asymptotic behav- 
ior, and study its properties under renormalization-group 
evolution, thereby obtaining strong constraints on model 
building. Using the results of this analysis, we propose a 
realistic model of the B-meson LCDA and use it to es- 
timate the important hadronic parameter Xb [2|, which 
enters in many analyses based on QCD factorization. 



II. MOMENT ANALYSIS 

The leading- twist, two-particle LCDA (fi® of the B- 
meson is defined in terms of the B-meson matrix ele- 
ment of a renormalized bilocal heavy-quark effective the- 
ory (HQET) operator relative to the matrix element of 
the corresponding local operator. The bilocal operator is 
made up of a soft spectator quark q s and a heavy quark 
h at light- like separation z, connected by a straight soft 
Wilson line S n (z, 0). Specifically, one defines 



[0\g s (z)S n {z,0)irh{0)\B{v)) 
(0\q s (0)iTh(0)\B(v)) 



(1) 



where r = v-z — ie. Our notation is such that z is propor- 
tional to a light-like vector n, v is the B-meson velocity, 
and for convenience we choose n ■ v — 1. The object T 
represents an arbitrary Dirac matrix chosen such that the 
operators have nonzero overlap with the B meson. The 
momentum-space LCDA is given by the Fourier trans- 
form 



1 

2^ 



(2) 



The analytic properties of the function 05 (r, /i) in the 
complex r plane imply that <f>^_ (uj, /x) = if u> < 0. 

We start by defining regularized moments of the B- 
meson LCDA as (for integer N > 0) 



M Ar (A uv ,/i) = 



duj uj n 4>+(uj, fx) . 



(3) 



A hard cutoff Atjv is imposed on the integral so as to 
avoid singularities from the region of large u> values, 
which are not regularized by renormalizing the bilocal 
operator in 0) [l| . The reason is that the position-space 
LCDA d)g(r, fi) and its derivatives are singular at t = 0. 
Only cut moments of the renormalized LCDA are UV fi- 
nite. For a sufficiently large value of Atjv the moments 
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Mjv(Auvj A*) can be expanded in a series of B-meson ma- 
trix elements of local HQET operators. The basic idea is 
the same as that used in previous work on cut moments 
of the B-meson shape function entering the analysis of in- 
clusive decays [3 HH ■ From the structure of the bilocal 
HQET operator in QJ and the Feynman rules of HQET 
it follows that the resulting local operators have Dirac 
structure 



(4) 



where the number of Dirac matrices inside the parenthe- 
sis is even if light quarks are treated as massless. By 
using the equations of motion iJj) q s = and iv ■ D h = 0, 
it is straightforward to find the corresponding operators 
of a given dimension D. For D — 3, the only possibility 
is the operator 



= q s t/i T h , 



(5) 



which appears in the denominator in (Q. For D = 4, 
there are naively four subleading operators with one 
derivative, namely 

Qia = Qs iv ■ D t/iT h , Qi c = q s i n ■ D i/iT h , 
Qib = q s in- Di/iTh, Qi d = q s ii/i p ijiT h . (6) 

However, the Wilson coefficients of the operators Q\ c and 
Qid are zero, because the residual momentum k of the 
external heavy-quark field only appears as v ■ k in HQET 
diagrams. Hence, these operators can be ignored. For 
D > 5, the situation becomes more complicated, since 
operators containing the gluon field need to be in- 
cluded. For our current analysis, we restrict the discus- 
sion to operators of dimension less than 5. 

The resulting expansion of the moments to subleading 
power in 1/ Ayv takes the form 



Mjv(A uv ,//) 



i—a.b 



A 



uv 



:0\Qn\B(v)) 
(0\Qa\5(v)) 



(7) 



where the ellipses denote terms of order (Aqcd/A 



uvj 



and higher. The short-distance coefficients Kn (AuVj A*) 
can be calculated using on-shell external quark states and 
employing partonic expressions for the LCDA and for 
the matrix elements of the local operators Q n to evalu- 
ate both sides of the matching relation Q . The relevant 
one-loop diagrams are shown in Figure^ Wave- function 
renormalization contributions cancel in the matching and 
thus can be omitted. We assign incoming residual mo- 
mentum k to the heavy quark and incoming momentum 
p to the light quark, subject to the on-shcll conditions 
v ■ k — and p 2 — 0. The Feynman amplitude is ex- 
panded to linear order in p before loop integrations are 
performed. This ensures that loop corrections to the 





1. 



FIG. 1: One-loop diagrams contributing to the partonic 
matrix elements of bilocal and local operators in HQET. A 
crossed circle denotes an operator insertion. Double lines rep- 
resent effective heavy-quark fields. 



matrix elements of the local operators vanish in dimen- 
sional regularization, because all integrals are scaleless. 
We thus obtain 

( | Qla |B(w))parton = V ■ p ( | Q \B(v)) 

parton j 

(0 1 Q lb |S(u)) parton =n-p(0\Q \B(v)) parton ■ (8) 

The result for the one-loop matrix element of the bilocal 
HQET operator is nontrivial. According to it pro- 
vides us with a partonic expression for the LCDA. After 
MS subtractions, we obtain at one-loop order 



A*) parton = 8(u)) \ l - 



CfOI s 7T 2 

12 



4-7T 

, Cfcxs 



in 

1 

w 



(9) 



2 - 



in 



-n ■ p 



-An ■ p 



CfOIs 



in 



1 + -r 



+ v ■ p 



CfCXs 



in 



In- 



+ (5n • p + 4v • p) 



where a s = a s (/j,) throughout, unless indicated oth- 
erwise. We have retained terms of linear order in p, 
which will be sufficient to extract the matching coeffi- 
cients Kq N ^ and K^P in J7J. The star distributions are 
generalized plus distributions defined as 

>] 



dui ■ 



F 



F In - , 
A* 



A 

dw F„ : 



In! 



= / dw ■ 



— ^ln^ 



Fo^A 
2 [i 



dw ■ 

Fo 
A 



F 



F^ln 



A 



(10) 



duiF^ 



In! 



= / duj ■ 
A 



Fn - w Fk 



ln 



A 



1 



2 /i 



where F(ui) is a smooth test function, and we use the 
short-hand notation F w = F(w) and F^ = F'{uj). 
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Given the results (JSJ) and it is straightforward to 
derive expressions for both sides of the matching relation 
(|7|) in the parton model, and to extract the desired ex- 
pressions for the Wilson coefficients. At one-loop order, 
we find 



K, 



(0) 



Ait 



-2 In 



A 



uv 



21n^V-^) 
H 12 J 



= (-4). 



K 



(0) 
lb 



4-7T 
47T 



4 In- 



A 



UV 



(11) 



for the zeroth moment, and 



K. 



(i) 



CfOL s 

Ait 



-AW- 



A 



uv 



A* 

Auv 



(12) 



"lb 1 



4-7T 



-2 In' 



A 



uv 



/' 



/i 12 



for the first moment. 

We have repeated the entire calculation outlined above 
in a different regularization scheme, where the depen- 
dence of the Feynman amplitudes on the component n ■ p 
of the light-quark momentum is kept exactly, whereas 
we linearize in the remaining components of p. In this 
scheme the loop corrections to the matrix elements of the 
local operators in |JSJ no longer vanish, and the result for 
the LCDA is far more complicated than that displayed in 
©. Nevertheless, we obtain the same expressions for the 

Wilson coefficients and K^P as given above. This 
is a highly nontrivial check, which gives us confidence in 
the correctness of our results. 

The Wilson coefficients describe the short-distance 
physics associated with the large cutoff scale Auv, an d 
hence it was legitimate to obtain them using a partonic 
calculation. Long-distance effects, on the other hand, 
reside in the hadronic matrix elements of the local oper- 
ators Q n , which cannot be calculated reliably using per- 
turbation theory. However, these matrix elements are 
constrained by heavy-quark symmetry and can be pa- 
rameterized in terms of universal form factors |26j . The 
results are particularly simple in the case of the operators 
Qu. Using relations derived in [2?]], we find that 



(0\Qia\B(v)) _ t (0\Qi b \B(v)) _AA 

(0|Qo|B(«)> ' (Q\Q Q \B(v)) 3 ' 1 ) 



where the quantity A = tub — rrib is the only hadronic 
parameter needed at this order. The first-order power 
corrections to the moments can now be expressed in 
terms of A. At one-loop order, and to subleading order 



in the power expansion in 1/A 



uv, 



M - 1 + 



Cfq>s 



An 

16A C F a s 



-2 In 



A 



uv 



M, = A 



3A UV 
Cfoi s 



Ait 



In 



Auv 



2 In 



- 1 



the results are 
Auv 1 



12 



4A 



uv 



1 



AlT 

C F a s 

47T 



-41n- 



A 



uv 



-21n 



Auv 



(14) 



iln- 



A 



uv 



7 

4 



r 2 \ 1 



12 



A 4 H 

These are our final expressions for the first two mo- 
ments of the renormalized I?-meson LCDA. As long as 
Auv 3* Aqcd, they are model-independent predictions 
of QCD, valid up to higher-order terms in a s and 1/Auv- 
The fixed-order perturbative expressions derived here are 
applicable if the two scales Auv and \i are of the same 
order, so that the logarithms in the matching coefficients 
are not parametrically large. 

Taking the derivative of the zeroth moment Mq in 
()14J) with respect to the cutoff, we can obtain a model- 
independent description of the asymptotic behavior of 
the B- meson LCDA El, i.e. 



B dM (A uv ,/x) 



dAuv 

At one-loop order, the result reads 



(15) 



Auv-w 



Cfoc s 



3ui 



ln^ 

H 



This relation holds for u> 3> Aqcd, 



(16) 

up to power cor- 



rections of order Aq CD /oj 3 . We observe that the radi- 
ation tail of the i3-meson LCDA becomes negative at 
lu w \/e/j, for a sufficiently large value of (i. This model- 
independent prediction for the asymptotic behavior of 
agrees qualitatively with the findings of the 
QCD sum-rule analysis in |23| . 



III. ELIMINATION OF THE POLE MASS 



Our calculations so far have been performed in the on- 
shell (pole) scheme, where A = raj — m^ olc is defined in 
terms of the 6-quark pole mass. However, it is well known 
that the pole mass suffers from infrared renormalon am- 
biguities [23, H3- Hence, it is desirable to eliminate 
the pole-scheme parameter A in favor of a new, short- 
distance parameter Ars defined in some renormalization 
scheme. For our purposes it is most convenient to employ 
a so-called "low-scale subtracted" heavy-quark mass de- 
fined with the help of a hard subtraction scale [i t . Exam- 
ples are the "kinetic mass" 1301 . the "potential-subtracted 
mass" 1^3, the "IS mass" |32j. and the "shape-function 
mass" [23,133]. Using the last definition as an example, 
we would use the relation 

A = A SF ^ f ,ri+^ f ^( 8 ln^-4)+... (17) 



4k 



/' 



4 



to eliminate the pole-scheme parameter A in the moment 
relations l|14fl. identifying the subtraction scale /if with 
the cutoff Auv- As always, a s = a s (fi). 

Alternatively, the moment relations themselves can be 
used to define a new subtraction scheme. Guided by the 
tree-level relations Mi = 4A/3 and Mo = 1, we are led to 
define a running parameter (the subscript "DA" stands 
for "distribution amplitude") 



AdaO/,^) 



3Afi(/x/,M) 
4M (m/,m) 



(18) 



to all orders in perturbation theory. From l|14|) , it follows 
that 



A = Ada(M/»m) 

CFOis 



47T 





^6 In 




-i)] 


An 








!)< 






V n 









(19) 



By taking the ratio of M\ and Mo in Ijl8(l the double- 
logarithmic radiative corrections are eliminated. Like the 
other short-distance mass definitions mentioned above, 
the parameter Ada can be regarded as a "physical" quan- 
tity in the sense that it is free of renormalon ambiguities. 
Perturbative relations can be used to transform from our 
new scheme to any other mass-definition scheme. For 
example, from l|17|) and l|19l) it follows that at one-loop 
order the parameter Ada is related to the parameter Asf 
in the shape-function scheme through 



Ada(m/,m) 



: A S f(m*,M*) 
C F a s 



-in 



1 + 



3 In 



4-7T 
El _ 

A 4 



6 In 



7 

A* 4 



(20) 



A rather precise value for Asf has been extracted from 
moment analyses of various spectra in the inclusive de- 
cays B — ► X s ~f and B — > X u lv, yielding Asf(m*, — 
(0.65 ± 0.06) GeV at /i* = 1.5 GeV (and at leading or- 
der in 1/mt) [33L l34j. This value will be used as an in- 
put when we compute the running parameter Ada (At/; A*) 
from the above relation. 



IV. RENORMALIZATION-GROUP 
EVOLUTION 

In Section^ we have derived model- independent pre- 
dictions for moments of the B-meson LCDA and for 
its asymptotic behavior for large w. The renormaliza- 
tion group can be used to obtain a model-independent 
description of how 0^(w, /i) changes under variation of 
the scale /i. The integro-differential evolution equation 
obeyed by the LCDA was derived in [35|, where an an- 
alytic solution was presented in the form of a double in- 
tegral. One finds that the distribution amplitude at a 



scale n can be expressed in terms of that at a lower scale 
A*o < M by 

1 [°° 



2tt 



where 



<Po(t) 



' du>' , R , , , / lo' 



(22) 



denotes the Fourier transform with respect to In to of the 
function 0^(o;,/Xo) at the initial scale jio- At leading 
order in perturbation theory, the kernel / takes the form 



/(a,,/x,/x ,i*) = e^°) (—) 
VMo/ 



it+g 



-27B9 



r(l-it-g)T(l + it) 

T(l + it + g)T(l- it) ' 1 J 



where 



and 







f da 


r C u SP (a) y 


J /?(«) 


a s («j) 





da 



7T + 7(a) 



(24) 



5 = .9(M,Mo) = 



(/ r cuS p(a) _ 2C F a s (^ ) 



a»(A»o) 



/3(a) 



/3 a s (/i) 



(25) 

In these expressions /3 = da, s /d\nu is the /?- function, 
and r cusp = Cpa s /TT + . . . , 7 = —Cpa s /2ir + . . . are 
anomalous dimensions. The perturbative expansion of 
V(n, /j,o) at next-to-leading order can be found in |3r| . 

Here we take a step further and simplify the solution 
obtained in |35| by performing the integration over t in 
(|21|l analytically. Substituting the expression for / from 
(|23|l . we observe that the integrand has poles situated on 
the imaginary axis in the complex t plane. The poles on 
the negative imaginary axis are located at t = —i(n — g) 
with n > 1 an integer (we assume < g < 1, which is 
satisfied for all reasonable values of scales), while those 
on the positive imaginary axis are located at t = in with 
n > 1 an integer. Using the theorem of residues, we 
obtain 



r(sO Jo ^ 



xl^l ^ 2Fl (l- 9 ,2- 9 ;2;^< 



UO J UJ> 



0J > 



(26) 



where uj < = min(aj,a/) and lo > — max(cj, to'). The hy- 
pergeometric function 2-Fi(a> b; c; 2) has the series expan- 
sion 



2 Fi(a,b; c;z) 



n=0 



(a)n (b)n Z n 

(c)„ n\ 



(27) 
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with (a)„ = L(a + n)/T{a). In the limit /i — > /iq we 
have V((j,,[io) -> 0, .9 -> 0, and 2 i<i(l - 5, 2 - 3; 2; x) -> 
(1 — x) 29_1 . Then the right-hand side in (|26|) reduces to 
the left-hand one. 

Equation H2(jf> provides the most compact expression 
possible for calculating the evolution of the LCDA under 
changes of the renormalization scale. It is tempting to 
conjecture that this is the exact solution to the evolution 
equation for the LCDA, valid to all orders in perturbation 
theory. An analogous statement is indeed true for the B- 
meson shape function |3?| |. In the present case, to prove 
this assertion one would need to show that the exact 
evolution equation for the LCDA is given by 



-jt 1" r cus P (a s ) In — + j(a s ) 

a In /i u 



r. 



j(a s ) / duj 

/o w > 



(28) 



TABLE I: Parameters of the model function (1291 for different 
values of the renormalization scale 



H [GeV] 


Ada [GeV] 


U)t [GeV] 


TV 


lj [GeV] 


1.0 


0.519 


2.33 


0.963 


0.438 


1.5 


0.635 


3.35 


0.974 


0.509 


2.0 


0.709 


4.32 


0.978 


0.557 


2.5 


0.770 


5.26 


0.981 


0.596 



substituting the model function l(2*9"|l into iJSJ, neglecting 



exponentially small terms 



-Auv/^o 



All remaining 



terms involving the cutoff Atjv are reproduced by con- 
struction, so that the results for N and ujq are indepen- 
dent of the cutoff, as they must be. At first order in a s , 
we obtain 



where L cusp is the universal cusp anomalous dimension 
of Wilson loops with light-like segments |3§, |3!| , and 7 
is some other anomalous dimension. In 35], the above 
relation was confirmed at one-loop order. 



V. PHENOMENOLOGICAL MODEL 

The model- independent properties of the J5-meson dis- 
tribution amplitude derived in this work provide use- 
ful constraints on model building. In this section we 
suggest a realistic form for (j)r_(ui, (i), which satisfies 
these constraints. For phenomenological purposes such 
a model is needed at a renormalization scale of order 
A* ~ \/ m b^QCD, as this is the characteristic "hard- 
collinear" scale for hard spectator scattering in exclu- 
sive B decays 0, |3^ . Our model consists of the two- 
component ansatz 

<fi(w, n) = N ^ e-W-o + 9{w - uh) — 



2 fi 



4A 



DA 



3lJ 



2 - In - 



(29) 



where Ada = Ada (/■*>/■*) i s defined in our new scheme 
(|19|l . and we set /if = /i for simplicity. The first term 
on the right-hand side is based on the exponential form 
proposed in p| , while the second piece is a radiation tail 
added so as to ensure the correct asymptotic behavior as 
shown in H6|) . The tail is "glued" onto the exponential 
at a position ui t chosen such that the resulting function 
is continuous. This yields 



y/efl ( 1 



2A 



DA 



14A DA 



\fejji 3e// 2 



(30) 



The normalization constant N and the parameter loq can 
be fixed by matching the expressions for the first two mo- 
ments in (|14f> with the corresponding results obtained by 



N = l 



= 1 



| C F a s 




21n 2 ^ 


47T 




A* 


16A DA 


In 




3u> t 




M / 


! Cfois 


(5 


TT 2 


An 




12 



UJf IT 

21n— - — 

fi 12 



and 



2A 



w = 



DA 



C F a s 



8A 



DA 



(31) 



Cfq s 

47T 

2A DA 



4tt 



L0 t 2 In 



61n^ 
A* 



16A 



DA 



3w t 
3/^ 



ln^-1 



(32) 



47T 



C F a s 

47T 



The expanded expressions for u>t, N, and loq are given 
only for the purpose of illustration. The exact expres- 
sions will be used in our numerical analysis. 

The model ansatz (|29|) has the attractive feature 
that it is to a good approximation invariant under 
renormalization-group evolution. Table [I] collects the pa- 
rameters entering this function for different values of /1, 
obtained using the central value Asf(m*j A 4 *) = 0-65 GeV 
in (I2U|I . For /i = 1 GeV and 2.5 GeV the correspond- 
ing functions are shown in Figure [21 For comparison, 
we also show the result at fx = 2.5 GeV obtained by ap- 
plying the evolution formula (|2bl to the model function 
at fi = 1 GeV. Both curves are very similar, indicating 
that the functional form (|29l) is approximately preserved 
under evolution. 

We have mentioned earlier that a QCD sum-rule anal- 
ysis of the B-meson LCDA at next-leading order in a s 
performed by Braun et al. |23| has exhibited an asymp- 
totic behavior similar to that of our perturbative QCD 
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FIG. 2: Model ansatz for the B-meson LCD A at fx = 1 GeV 
(narrow solid curve) and 2.5 GeV (wide solid curve). The 
dashed curve shows the result at 2.5 GeV obtained by evolving 
the distribution amplitude from 1 to 2.5 GeV. 




FIG. 3: Two different models for the B-meson LCDA at 
/i = 1 GeV, constrained to have the same normalization and 
first moment. The solid curve corresponds to (1291 . the dashed 
one to 1331 1. 



analysis. These authors have proposed the model form 

1 2(a B - 1) 



k 

k 2 + l 



1 



Ink 



(33) 

at fx = 1 GeV, where k = w/1 GeV. The two parame- 
ters entering this functions are defined in terms of the 
integrals 



A. 1 



a B \ B l = - / 
Jo 



aui — 



In- 



(34) 



The parameter ranges obtained from the sum-rule anal- 
ysis arc A^ 1 = (2.15 ± 0.50) GeV" 1 and a B = 1.4 ± 0.4 
at \i = 1 GeV. On the other hand, if we require that 
the function obey the moment constraints l|14|) at a 
large value of the cutoff, say Atjv = 3 GeV, then we find 
Xg 1 = (1.79 ±0.06) GeV" 1 and a B = 1.57 ±0.27. These 




1 2 3 4 5 6 

A UV [GeV] 

FIG. 4: Comparison of model results (black) and OPE predic- 
tions (gray) for the first two moments of the LCDA, evaluated 
at /i = 1 GeV and for different values of the cutoff. The solid 
black curves are obtained in our model 12911 . the dashed ones 
in the model (1331 of |2^ . 



values are consistent with the findings of |2J| • It is inter- 
esting that, once the moment constraints are imposed, 
the two models in (|29(l and l|33|l are nearly indistinguish- 
able, in spite of the rather different functional forms (ex- 
ponential vs. power-like fall-off). This fact is illustrated 
in Figure |3I 

We are now in a position to investigate how moments 
of the LCDA computed using the model functions l|29(l 
and (|33|l compare with the model-independent predic- 
tions (fTlj) of the OPE, which are valid for Atjv ^ Aqcd- 
In Figure 01 we show in black the model results for the 
moments Mq and Mi at /x = 1 GeV as a function of 
the cutoff Aijv- F° r comparison, the gray curves show 
the predictions of the OPE. We observe that our model 
curves quickly converge toward the OPE predictions for 
Atjv > 2.5 GeV. For large cutoff values the agreement is 
perfect, since by construction our function has the cor- 
rect asymptotic behavior. The model of Braun et al. 
agrees qualitatively with the OPE for large Atjv, but ex- 
act agreement can only be enforced at a single value of 
the cutoff (3 GeV in our case) . Note that for small val- 
ues of Atjv there are significant deviations between the 
OPE predictions and the model results. This is expected, 
given that the OPE is only valid for Atjv 3> Aqcd- For 
Atjv = 2 GeV, for example, we expect unknown cor- 
rections of order (A/Atjv) 2 ~ 0.1 to M , and of order 
A 2 /A uv ~ 0.2 GeV to M x . This is consistent with the 
deviations seen in the figure. 
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TABLE II: Inverse moments A s * and erg calculated using 
the model function 12911 



H [GeV] 


As 1 [GeV- 1 ] a B 


1.0 


2.09 ±0.24 1.61 ±0.09 


1.5 


1.86 ±0.17 1.79 ±0.08 


2.0 


1.72 ±0.14 1.95 ±0.07 


2.5 


1.62 ±0.12 2.09 ±0.07 



VI. ESTIMATES FOR INVERSE MOMENTS 

The "inverse moments" defined in l|34l) play an im- 
portant role in the analysis of many exclusive B-meson 
decays. They control the strength of the leading-power 
spectator interactions in leptonic decays such as B — > 
jlv, semileptonic decays such as B — > ttIu, and hadronic 
decays such as B — > ixix. The quantity as enters these 
analyses as soon as one goes beyond the tree approx- 
imation. Given that we have constructed highly con- 
strained models for the distribution amplitude which sat- 
isfy the QCD predictions for moments and have the cor- 
rect asymptotic behavior, it is interesting to ask what 
estimates we can obtain for the parameters Xb and erg. 

In Table |n] we collect the results for the two inverse 
moments obtained using the model ansatz (|29() . The er- 
ror bars reflect the variation of the results with the input 
parameter Asf = (0.65 ± 0.06) GeV. In addition, there 
are other theoretical uncertainties related to the neglect 
of higher-order terms in the OPE and, more importantly, 
to nonperturbative hadronic uncertainties in the precise 
shape of the LCDA for small values of ui. For instance, 
comparing the results in the table with those obtained 
using the model (|33[1 at ju = 1 GeV, we observe shifts 
in A^, 1 and as by 0.3 GeV^ 1 and 0.04, respectively. We 
believe that the true theoretical uncertainties are about 
twice as large as the errors shown in the table. A graph- 
ical representation of the results is shown in Figure 
where the light gray bands are an estimate of the total 
theoretical uncertainty. 

Our findings are in good agreement with the QCD sum- 
rule estimates at next-to-leading order in a s obtained by 
Braun et al. |2^|, indicated by the data points in the 
figure. We may also compare with earlier estimates of 
Xg 1 derived from lowest-order QCD sum rules, where 
the scale dependence is not controlled. Grozin et al. [| 
found Xg 1 = 3/(2A) « 2.2 GeV" 1 (for a typical value 
m b « 4.6 GeV), while Ball et al. obtained A^ 1 « 
1.7 GeV. Both are consistent with our findings. 



VII. CONCLUSIONS 

Using rigorous methods based on the operator prod- 
uct expansion, we have studied some model-independent 




1 

1 1.5 2 2.5 3 3.5 4 

fi [GeV] 

3 



2.5 

cq 2 
b 

1.5 




1 U 1 

1 1.5 2 2.5 3 3.5 4 

fi [GeV] 

FIG. 5: Model estimates of the inverse moments A^ 1 and 
aB for different values of the renormalization scale. The dark 
bands reflect the uncertainty in the value of A, whereas the 
light bands represent an estimate of the total theoretical error. 
The data points show the results obtained from the QCD sum- 
rule analysis of [2^1 . 



properties of the £?-meson light-cone distribution am- 
plitude <fiB(ui,fi). We have derived explicit expressions 
for the first two moments of the distribution amplitude 
as a function of the renormalization scale /i and a hard 
Wilsonian cutoff Atjv applied to integrals over u>. The 
ratio Mi/Mo of the first two moments can be used to 
define a physical subtraction scheme for the parameter 
A = tub — mb of heavy-quark effective theory. This 
links the only nonperturbative hadronic parameter en- 
tering the moment predictions at next-to-leading power 
in 1/Atjv in a calculable way to the 6-quark mass. From 
the cutoff dependence of the moment Mo we have derived 
an analytic expression for the asymptotic behavior of the 
distribution amplitude for large uj 3> Aqcd , valid at first 
order in a s and at next-to- leading order in 1/lu. Finally, 
we have presented a new, compact evolution formula that 
expresses the distribution amplitude at some scale /i in 
terms the function </>+(w, fio) at a lower scale /io- 

Based on our analysis we have proposed a realistic 
model of the B-meson distribution amplitude, which is 
consistent with the moment relations. With the help of 
this function we have obtained estimates for the inverse- 
moment parameters A b and aB , which play an important 
role in many phenomenological applications of the QCD 
factorization approach to exclusive B decays. We find 
A^ 1 = (2.1±0.5) GeV- 1 and a B = 1.6±0.2 at fi = 1 GeV 
with conservative errors. 
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We hope that our analysis will not only supply a guide- 
line for understanding the i?-meson distribution am- 
plitude without relying on a specific model, but also 
open a new strategy for further, more detailed studies 
of </)+(uj,(j,) using a systematic short-distance approach. 
Ultimately, this may help to reduce the theoretical un- 
certainties in predictions for exclusive i?-meson decays. 
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